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Detection of visible light from the darkest world 
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ABSTRACT 

We present the detection of visible light from the planet TrES-2b, the darkest exoplanet 
currently known. By analysis of the orbital photometry from publicly available Kepler 
data (0.4-0.9 ^m), we determine a day-night contrast amplitude of 6.5 ± 1.9 ppm, 
constituting the lowest amplitude orbital phase variation discovered. The signal is 
detected to 3.7cr confidence and persists in six different methods of modelling the data 
and thus appears robust. In contrast, we are unable to detect ellipsoidal variations or 
beaming effects, but we do provide confidence intervals for these terms. If the day- 
night contrast is interpreted as being due to scattering, it corresponds to a geometric 
albedo of A g = 0.0253±0.0072. However, our models indicate that there is a significant 
emission component to day-side brightness, and the true albedo is even lower (< 1%). 
By combining our measurement with Spitzer and ground-based data, we show that 
a model with moderate redistribution (P„ ~ 0.3) and moderate extra optical opacity 
(«/ ~ 0.3 — 0.4) provide a compatible explanation to the data. 
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1 INTRODUCTION 

Orbital photometric phase variations have long been used in 
the s tudy and characterisation of eclipsing binaries l|Wilsonl 
1 1994 ). where the large masses and small orbital radii re- 
sult in phase variations at the magnitude to millimagnitude 
level. The three dominant components of these variations 
are i) ellipsoidal variations, due to the non-sph erical nature 
of a s tar caused by gravitational distortion (e.g. IWelsh et alj 
|2010| ) ii) relativistic beaming, due to th e radial motion of the 
star shifting the stellar spectrum (e.g. iMaxted et alJl2~00Ch 
iii) reflected/emitted light, wh ich varies depen ding on what 
phase of a body is visible (e.g. iFor et al.ll2010l ). 

The visible bandpass orbital phase variations of a star 
due to a hot- Jupiter companion are much smaller - around 
the part-per-million (ppm) level - and thus have eluded de- 
tection until relatively recently. The high pre cision space- 
based photometry of C0R0T (0.56-0 .71 (im) l|Baglin et al.l 
120091 ) and Kepler (|Basri et al.l 120051 ) have opened up this 
exciting new way of studying exoplanets for first time, with 
several detections recently reported in the literature: 



■ CoRoT-lb (|Snellen et al.l|2009| ); reflected/emitted light 
amplitude 126 ± 36 ppm 

■ HAT-P-7b (|Welsh et all l201Ch : ellipsoidal amplitude 
37 ppm, reflecte d /emitted light amplitu de 63.7 ppm 

■ CoRoT-3b (|Mazeh fc Faiglerl l20ld ); ellipsoidal ampli- 
tude (59 ± 9) p pm, beaming amplit ude (27 ± 9) ppm 

■ Kepler- 7b l|Demorv et al.ll201ll ); reflected/emitted light 
amplitude (42 ± 4) ppm 



In this letter, we investigate the hot- J upiter orbiting the 
GOV star TrES-2 i|Q'Donovan et al.ll2006l ). where we detect 
a reflected/emitted light amplitude of (6.5 ± 1.9) ppm to a 
confidence of 3.7cr, or 99.98%. We also measure the ellip- 
soidal variation and relativistic beaming amplitudes to be 
(1.5 ± 0.9) ppm and (0.2 ± 0.9) ppm respectively, which are 
broadly consistent with theoretical expectation. 

If our detected signal is interpreted as being purely 
due to scattering, then the corresponding geometric albedo 
would be A g = 0.0253 ± .0072 (using system pa rameters 
from Table 2, column 2 of iKipping fc Bakosl [201ll (KB11), 
as will be done throughout this work), meaning that just 
four months of Kepler's exquisite photometry has detected 
light from the darkest exoplanet yet found. Extrapolating to 
a 6 year baseline, one can expect to detect albedos 0.1 (to 
3a confidence) at similar orbital radii down to Rp ~ 3.0 R®. 
This clearly highlights the extraordinary potential which 
would be granted by an extended mission for Kepler. 

2 OBSERVATIONS & ANALYSIS 
2.1 Data Acquisition 

We make use of "Data Release 3" (DR3) from the Kepler- 
Mission, which consists of quarters 0, 1 and 2 (Q0, Ql & 
Q2). Full details on the data processing pipeline can be 
found in the DR3 handbook. The data includes the use of 
BJD (Barycentric Julian Date) time stamps for each flux 
measurement, which is crucial for time sensitive measure- 
ments. All data used are publicly available via MAST. 
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We use the "raw" (labelled as "AP_RAW_FLUX" in 
the header) short-cadence (SC) data processed by the DR3 
pipeline and a detailed description can be found in the ac- 
companying release notes. The "raw" data has been subject 
to PA (Photometric Analysis), which includes cleaning of 
cosmic ray hits, Argabrightenings, removal of background 
flux, aperture photometry and computation of centroid po- 
sitions. It does not include PDC (Pre-search Data Condi- 
tioning) algorithm developed by the DAWG (Data Analysis 
Working Group). As detailed in DR3, this data is not rec- 
ommended for scientific use, owing to, in part, the potential 
for under/over- fitting of the systematic effects. 

2.2 Cleaning of the Data 

The raw data exhibit numerous systematic artifacts, includ- 
ing pointing tweaks (jumps in the photometry), safe mode 
recoveries (exponential decays) and focus drifts (long-term 
trends). The first effect may be corrected by applying an off- 
set surrounding the jump, computed using a 30-point inter- 
polative function either side. Due to the exponential nature 
of the second effect, we chose to exclude the affected data 
rather than attempt to correct it. The third effect may be 
corrected for using a detrending technique. 

For this latter effect, we use the cosine filter 
utilised to detect elli psoidal variations in CoRoT data by 
iMazeh fc Faigleil (|2010h . The technique acts as a high-pass 
filter allowing any frequencies of the orbital period or higher 
through and all other long-term trends are removed. Thus, 
we protect any physical flux variations on the time scale 
of interest. We applied the filter independently to Q0+Q1 
data and then Q2 data. This is because the Kepler space- 
craft was rotated in the intervening time and so the long- 
term trend will not be continuous over this boundary. After 
removing 3a outliers with a running 20-point median and 
transits using the ephemeris of KB11, we applied the filter, 
with the resulting fitted trends shown in Fig [1] Our final 
cleaned data consists of 154,832 SC measurements with a 
mean SNR~ 4408. 

2.3 Three Models 

We first define our null hypothesis, Mi, where we employ a 
flat line model across the entire time series, described by a 
constant ao. For a physical description of the orbital phase 
variations, we first tried the same model as that us ed by 
ISirko fc Paczvnskil l|2003h and lMazeh fc Faigleri l|2010r i. This 
simple model is sufficient for cases where one is dealing with 
low signal-to-noise and reproduces the broad physical fea- 
tures. The model, M2, is given by 

M.2(4>) = ao + 0.5ai c cos(<t5) + ai s sin(<7j) 

+ a 2c cos(20) + a 2s sin(2c» , (1) 

where (j> is the orbital phase (defined as being at the 
time of transit minimum) and a,i are coefficients related to 
the physical model, ao is simply a constant to remove any 
DC (direct-current) component in the data. a± c corresponds 
to the reflection/emission effect and is expected to be have a 
negative amplitude. a\ s corresponds to the relativistic beam- 
ing effect and is expected to be positive. a 2c corresponds to 
the ellipsoidal variations and should be negative. a 2s is a 
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Figure 1. "Raw" (PA output) flux from DR3 of the Kepler 
pipeline for Q0&Q1 (upper panel) and Q2 (lower panel) of the 
star TrES-2. Overlaid is our model for the long-term trend, com- 
puted using a discrete cosine transform for each data set. Outliers 
and discontinuous systematic effects have been excluded. 



dummy term which should be zero and ensures the ellip- 
soidal variation is detected with the correct phase. 

We also tried a third model, where the ai c term is 
replaced by the reflection caused by a Lambertian sphere: 



0.5ai c cos < 



die 



+ (n — |0|) cos I 



(2) 



2.4 Three Data Modes 



In addition to three models, we have three data input modes. 
The first is simply corrected for detrending and nothing else, 
denoted T>i . The second mode renormalises each orbital pe- 
riod epoch. This renormalisation is done by computing the 
median of each epoch and dividing each segmented time 
series by this value and we denote this mode as T>2. Fi- 
nally, we tried allowing each orbital period epoch to have its 
own variable renormalisation parameter, which is simulta- 
neously fitted to the data along with the orbital phase curve 
model. This parameter is dubbed ao,j for the j th orbital pe- 
riod epoch. Denoting this data input mode as D3, the fits 
now include an additional 51 free parameters. 

The models are fitted to the unbinned data using a 
Markov Chain Monte Carlo algorithm described in KB11 
(method A) with 1.25 x 10 accepted trials burning out the 
first 25,000. In total, there are nine ways of combining the 
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three models with the three data modes. All nine models are 
fitted and results are given in Table [TJ with our preferred 
model description being M2, T>3 (since thermal emission is 
likely dominant over scattering, see |4]). 

3 RESULTS 

3.1 Orbital Photometry 

Table [T] presents the results of fitting the detrended Kepler 
photometry. Our models make no prior assumption on the 
sign or magnitude of the a,i coefficients. The orbital period 
and transit epoch are treated as Gaussian priors from the 
circular orbit results of KB11. 

When considering statistical significance, what one is 
really interested in is the confidence of detecting each phys- 
ical effect i.e. reflection/emission, ellipsoidal and beam- 
ing. For this reason, model comparison tools, such as the 
Bayesian Information Criterion (BIC) or an F-test are inap- 
propriate. This is because these methods evaluate the pref- 
erence of one hypothesis over another, where the two mod- 
els would be a null-hypothesis and a hypothesis including 
reflection/emission, ellipsoidal variations, beaming and the 
dummy term. Thus, any inference drawn from this would be 
for the entire model and not for each individual effect. In the 
analysis presented here, simple inspection of the posteriors 
from Fig[2]shows that only one effect is actually detected (re- 
flection/emission), but a model comparison method would 
evaluate the significance of all four physical effects (includ- 
ing the dummy term) versus no effect. 

A more useful statistical test would consider the signif- 
icance of each physical effect individually from a joint fit. 
A n excellent tool to thi s end is the odds ratio test discussed 
in iKipping et all (|201Cf ). If a parameter was equal to zero, 
we would expect 50% of the MCMC runs to give a positive 
value and 50% to give a negative value. Consider that some 
asymmetry exists and a fraction / of all MCMC trial were 
positive and 1 — / were negative. The reverse could also be 
true and so we define / such that / > 0.5 i.e. it represents 
the majority of the MCMC trials. The odds ratio of the 
asymmetric model over the 50:50 model is: 

For only two possible models, the probability of the 
asymmetric model being the correct one is P(asym) = 
1 — [1/(1 + O)]. We perform this test on each of the four 
parameters fitted for, ai c , 02e, ai s and 02 S . The associated 
results are visible in the top-left corners of each posterior 
shown in Fig [2] for our preferred model and data mode i.e. 
M2, 2?3. To summarise, only one parameter presents a sig- 
nificant detection - the reflection/emission effect. Here, we 
find die's posterior is sufficiently asymmetric to have a prob- 
ability of occurring by random chance of just 0.02%, which 
equates to 3.67a. We consider any signal detected above 3a 
confidence to merit the claim of a "detection" rather than 
a measurement and thus we find TrES-2b to be the darkest 
exoplanet from which visible light has been detected. 

As discussed in §2.31 we tried two different models for 
the reflection/emission effect; a simple sinusoid (M2) and 
the reflected light from a perfectly Lambertian sphere (M3). 
Between the two models, there is negligible difference in 



the goodness-of-fit, as seen in Table [T] for all three data 
modes. Including the Lambertian model takes some power 
away from the ellipsoidal variations though and thus the cur- 
rent data does not yield a preference between a Lambertian 
sphere model or stronger ellipsoidal variations. 

3.2 Occultation Measurement 

The duration of the transit, and thus occultation since 
TrES-2b maintains negligible eccentricity, is equal to 4624 ± 
42 seconds (defined as the time between when the planet's 
centre crosses the stellar limb to exiting under the same 
condition). In contrast, the orbital period of TrES-2b is 
2.470619 days. We therefore obtain ~46 times more integra- 
tion time of the orbit than the occultation event. This indi- 
cates that we should expect to be able to reach a sensitivity 
of \/46 times greater, purely from photon statistics. The un- 
certainty on our phase curve measurement is ±1.9 ppm. We 
therefore estimate that one should have an uncertainty on 
the occultation depth of ~ 13 ppm. If we assume the night- 
side has a negligible flux, then the depth of the occultation 
is expected to be 6.5 ppm (i.e. equal to the day-night con- 
trast), and this already suggests that the present publicly 
available Kepler photometry will be insufficient to detect 
the occultation. To test this hypothesis, we will here fit the 
occultation event including the Q0, Ql and Q2 data. 

To perform our fit, we use the same Gaussian priors on 
P and r as earlier. We also adopt priors for other important 
system parameters from KB11, such as b = 0.8408 ± 0.0050, 
p 2 = 1.643 ± 0.067% and fx = 4624 ± 42 seconds. We stress 
that these are all priors and not simply fixed parameters. 
We also make use of the priors on the ao,j coefficients from 
the M3, D3 fit. Data are trimmed to be within ±0.06 days 
of the expected time of occultation to prevent the phase 
curve polluting our signal, leaving us with 8457 SC mea- 
surements. Assuming a circular orbit, the data were fitted 
using an MCMC algorithm. 

The marginalised posterior of the occultation depth 
yields 5 occ = 16 1.14 ppm, which is clearly not a significant de- 
tection. The derived uncertainty of 13-14 ppm is very close to 
our estimation of ~ 13 ppm and thus supports our hypothe- 
sis that the current Kepler data are insufficient to detect the 
occultation of TrES-2b. We also note that the inclusion of 
the Q2 data does improve the constraints on the occultation 
event (KB11 found Socc = 2ll|| ppm using Q0 & Ql only). 

4 DISCUSSION 

Hot-Jupiters are generally expected to be dark. Significant 
absorption due to the broad wings of the sodium and potas- 
sium D lines is thou ght to dominate their visible spectra 
l|Sudarskv et al.ll2000l ). leading to low albedos of a few per- 
cent. Indeed, aside from the recent report of Kepler-7b's 
(38 ±12)% Kepler-band geometric albedo (KBU), searches 
for visible light fr om hot-Jupiters have generally revealed 
mere upper limits ([C ollier Cameron et al. 2002; Lcigt i"et al.l 
l2003l ; lRowe et al.ll2008l ; iBurrows et al.ll2008h . 

The 6.5±1.9ppm contrast (determined from our pre- 
ferred model M2, T>s) between the day-side and night-side 
photon flux that we measure for TrES-2b represents the 
most sensitive measurement yet of emergent radiation in 
the visible from a hot- Jupiter, and is a factor of ~20 and 
~6 dimmer than the corresponding differences for HAT-P- 
7b l|Welsh et alj|20ich and Kepler-7b. 
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Table 1. Results of three models with three data modes, giving nine sets of results. Emboldened row denotes our favoured solution. Results 
do not include the orbital period, P and transit epoch r, which are treated as Gaussian priors via P = 2.47061896 ± 0.00000022 days and 
t = 2454950.822014 ± 0.000027 BJDtdb • Quoted values are medians of marginalised posteriors with errors given by la quantiles. * = 
parameter was fixed. We do not show the aoj fitted terms, which are simply renormalisation constants and are available upon request. 
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Figure 2. Top Left: Final fit to the phased photometry. Points without errors are the 2000-point phase binned data. Points with 
errors are 5000-point phase binned data. Best-fit model M2 with data mode T>z shown in solid. Note that all fits were performed on the 
unbinned photometry. Top Right & Lower Panels: Marginalised posterior distributions for the same model of four fitted parameters. 
Unity minus the false-alarm-probability values are provided for each parameter, based upon an odds ratio test described in $3.1\ 



In order to interpret the visible flux, we use the plane- 
tary atmosphere modelling code C00LTLUSTY (|Hubenv et al.l 
120031 ). For simplicity, we adopt equilibrium chemistry with 
nearly Solar abundance of elements, although we leave tita- 
nium oxide and vanadium oxide (TiO and VO) out of the 
atmosphere model. These compounds could, if present in the 
upper atmosphere of a hot- Jupiter, strongly affect the atmo- 
sphere structure and the visible and near infrared spectra, 
by making the atmosphere more opaque in the visible and 
by leading to a thermal in version if the stellar irradiation ex- 
ceeds ~10 9 erg cm -2 s" 1 (|Hubenv et al.ll2003l ; lFortnev et all 
2008). We leave TiO and VO out of our calculations, how- 
ever, because of the difficulty of maintaining heavy, conden- 



sible species high in the atmospheres l|Spiegel et alj|200^ ). 
Instead, we use an ad hoc extra op acity source k! , as de- 
scribed in ISpiegel fc Burrows! (|2010l ). 

We calculate a grid of models with k' ranging (in 
cm 2 g _1 ) from to 0.6 in steps of 0.1 and with redistri- 
bution P n ranging from to 0.5 in steps of 0.1 (P„ repre- 
sents the fraction of incident irradiation that is transported 
to the nightside, which is assumed in our models to occur 
in a pressure range from 10 to lOOmbars). For each of these 
42 parameter combinations, we calculate a day-side model, 
a night-side model, and a model that has the same tempera- 
ture/pressure structure as the dayside but that has the star 



© 2011 RAS, MNRAS 000, [T]{5] 



The Darkest World 5 



turned off, so as to calculate the emitted (and not scattered) 
flux (thus also giving the scattered component). 

We draw several inferences from our models and the 
data. First, the nightside contributes negligible flux in 
the Kepler-band (always <12% of the dayside, and for 
most models significantly less than that), meaning that the 
6.5 ppm number represents essentially the entire day-side 
flux. 

Second, by also including the available infrared sec- 
ondary eclipse d ata on TrES-2b (|Q'Donovan et all l20ld ; 
ICroll et al.|[2010f ), we find that in our model set there must 
be some redistribution (but not too much) and there must be 
some extra absorber (but not too much). For each model, 
we compute a \ 2 value, including 6 data points: Kepler- 
band, _R"s-band, and the four Spitzer IRAC channels (3.6, 
4.5, 5.8 & 8.0/xm). Fig portrays the \ 2 values of our 
grid of models, with the colour ranges corresponding to 
the x 2 values bounding 68.3% of the integrated probabil- 
ity (la), 95.5% (2a), 99.7% (3<r) and 99.99% (4a). The 
models that best explain the available data correspond to 
ft' ~ 0.3 - 0.4 cm 2 g" 1 and P n ~ 0.3 (~30% of incident 
flux redistributed to the night). In particular, models with 
no extra absorber are completely inconsistent with obser- 
vations, even on the basis of the Kepler data alone. The 
upshot is that some extra opacity source appears to be re- 
quired to explain the emergent radiation from this extremely 
dark world. Owing to this optical opacity, our models that 
are consistent with the data have thermal inversio ns in their 
upp er atmosphere, as in | Spiegel fc Burrows! (|2010l ) . We note 
that lMadhusudhan fc Seaeeri (j2010f ) find that the IR data of 



TrES-2b may be explained by models both with and with- 
out thermal inversions; nevertheless, we believe that optical 
opacity sufficient to explain t he Kepler data i s likel y to heat 
the upper atmosphere, as per iHubenv et all (|2003h . 

Finally, by computing the scattered contribution to the 
total flux, we find that for all parameter combinations the 
scattered light contributes <10% of the Kepler-band flux, 
and for the best-fit models the scattered light is <1.5% of 
the total. TrES-2b, therefore, appears to have an extremely 
low geometric albedo (for all models, the geometric albedo 
is < f%, and for the best-fit models it is ~0.04%). Exact 
values for the amount of extra optical opacity, redistribu- 
tion and the albedo cannot be presently provided because 
inferences about them depend on unknown quantities such 
as the wavelength dependence of the extra opacity source 
and the altitude dependence of winds. 

ACKNOWLEDGMENTS 

We thank the Kepler Science Team, especially the DAWG, 
for making the data used here available. Thanks to A. Bur- 
rows, M. Nikku & the anonymous referee for helpful com- 
ments and I. Hubeny & A. Burrows for the development 
and continued maintenance of COOLTLUSTY and associated 
opacity database. DMK is supported by Smithsonian Instit. 
Restricted Endowment Funds. 

REFERENCES 

Baglin, A. et al. 2009, Transiting Planets, Proc. IAU 

Symp., 253, 71 
Basri, G., Borucki, W. J. & Koch, D. 2005, New Astronomy 

Rev., 49, 478 



0.6 












I 


0.5 






UU 

















j 


E 0.3 












1 


* 0.2 














0.1 












□ 





■■ 


■ 


■ 


■ 


□ 



3-4o 



1-2o 



o 

o 
o 

CL 

CD 
cn 
cn 

o 

Tl 



0.1 0.2 



0.3 0.4 
P 

n 



0.5 



Figure 3. Goodness- of -fit for a grid of atmosphere models. The 
models that are consistent with available Kepler-ftartd, Ks-band, 
and Spitzer IRAC data have moderate redistribution to the night 
side (P n ) and moderate extra optical opacity (k! ). Models with 
k 1 = can be ruled out on the basis of the Kepler data alone. 



Burrows, A., Ibgui, L., & Hubeny, I. 2008, ApJ, 682, 1277 

Demory, B.-O. et al. 2011, ApJL, accepted 

Collier Cameron, A., Home, K., Penny, A. & Leigh, C. 

2002, MNRAS, 330, 187 
Croll, B., Albert, L., Lafreniere, D., Jayawardhana, R., & 

Fortney, J. J. 2010, ApJ, 717, 1084 
For, B.-Q., et al. 2010, ApJ, 708, 253 

Fortney, J. J., Lodders, K., Marley, M. S., & Freedman, 

R. S. 2008, ApJ, 678, 1419 
Hubeny, I., Burrows, A., & Sudarsky, D. 2003, ApJ, 594, 

1011 

Kipping, D. M. et al. 2010, ApJ, 725, 2017 

Kipping, D. M. & Bakos, G. A. 2011, ApJ, 733, 36 (KB11) 

Leigh, C, Collier Cameron, A., Home, K., Penny, A. & 

James, D. 2003, MNRAS, 344, 1271 
Madhusudhan, N., & Seager, S. 2010, ApJ, 725, 261 
Maxted, P. F. L., Marsh, T. R. & North, R. C. 2000, MN- 
RAS, 317, L41 
Mazeh, T. & Faigler, S. 2010, A&A, 521, 59 
O'Donovan, F. T. et al. 2006, ApJ, 651, L61 
O'Donovan, F. T., Charbonneau, D., Harrington, J., Mad- 
husudhan, N., Seager, S., Deming, D. & Knutson, H. A. 
2010, ApJ, 710, 1551 
Rowe, J. F. et al. 2008, ApJ, 689, 1345 
Sirko, E. & Paczynski, B. 2003, ApJ, 592, 1217 
Snellen, I. A. G., de Mooij, E. J. W. & Albrecht, S. 2009, 

Nature, 459, 543 
Spiegel, D. S., Silverio, K. & Burrows, A. 2009, ApJ, 699, 
1487 

Spiegel, D. S. & Burrows, A. 2010, ApJ, 722, 871 
Sudarsky, D., Burrows, A., & Pinto, P. 2000, ApJ, 538, 885 
Welsh, W. F., Orosz, J. A., Seager, S., Fortney, J. J., Jenk- 
ins, J., Rowe, J. F., Koch, D. & Borucki, W. J. 2010, ApJ, 
713, 145 

Wilson, R. E. 1994, PASP, 106, 921 



© 2011 RAS, MNRAS 000, [TH5] 



